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1 INTRODUCTION 



ABSTRACT 

We study the variation of the dark matter mass fraction of elliptical galaxies as a 
function of their luminosity, stellar mass, and size using a sample of 29,469 elliptical 
galaxies culled from the Sloan Digital Sky Survey. We model ellipticals as a stellar 
Hernquist profile embedded in an adiabatically compressed dark matter halo. This 
model allows us to estimate a dynamical mass (Mdynm) at the half-light radius from the 
velocity di spersion of the spectra, a nd to compare these to the stellar mass estimates 
(M*) from lKauffmann et alJ (|2003f) . We find that M*/L is independent of luminosity, 
while Mdynm/L increases with luminosity, implying that the dark matter fraction 
increases with luminosity. We also observe that at a fixed luminosity or stellar mass, 
the dark matter fraction increases with increasing galaxy size or, equivalently, increases 
with decreasing surface brightness: high surface brightness galaxies show almost no 
evidence for dark matter, while in low surface brightness galaxies, the dark matter 
exceeds the stellar mass at the half light radius. We relate this to the fundamental plane 
of elliptical galaxies, suggesting that the tilt of this plane from simple virial predictions 
is due to the dark matter in galaxies. We find that a simple model where galaxies are 
embedded in dark matter halos and have a star formation efficiency independent of 
their surface brightness explains these trends. We estimate the virial mass of ellipticals 
as being approximately 7-30 times their stellar mass, with the lower limit suggesting 
almost all of the gas within the virial radius is converted into stars. 

regularities is still an open problem; they provide a strong 
constraint and challenge to theories of galaxy formation. 



The zoo of galaxies is an eclectic one, with large variations in 
morphologies, colours and spectra. This is perhaps not very 
surprising, given the varied environments in which galaxies 
reside and the different processes underlying their forma- 
tion. However, a number of regular trends are also known; 
for instance, the rotation velocities of spiral galaxies and 
the central velocity dispersions of e lliptical galaxies corre- 
late strongly with luminosit y - the iTullv fc Fisherl Jl977ft 
and iFaber fc Jackson] jl976Tl relations. In addition, ellipti- 
cal galaxies appear to reside in a thin plane in the space 
of their luminosity (L), radius (R), and central velocity dis- 
persion (a), the "fundam ental plane" dDiorgovski fc Davisl 
Il987t iDressler et alJll987h . of which the Faber- Jackson re- 
lation is a projection. Understanding the source of these 



Elliptical galaxies are probably the simplest galaxies to 
model. These are pressure supported systems, and are re- 
markably uniform in their properties. Their intens ity pro- 
files are well described by a Ide Vaucouleurs! lll948t) profile, 
I oc exp(-(6>/6l s ) 1/4 ), they have very simple spectral energy 
distributions and they have very uniform photometric prop- 
erties. However, very little has been known about the mat- 
ter content of elliptical galaxies until recently, principally 
due to the lack of dynamical tracers like the HI rotation 
curves for spiral galaxies. More recently however, there have 
been a number of surveys that have measured the dynami- 
cal struct ure of samples of elliptica l galaxies using slit spec- 
troscopy dKronawitter et alj|2000 |), and with the advent of 
integral field spectroscopy ide Zeeuw et aUl2002lh the num- 
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ber of well studied elliptical galaxies will only increase. In ad- 
dition , strong lensing measurements fee. iKoopmans fc Treul 
2003) have added independent measures of the masses of the 
halos of ellipticals, especially at radii outside the realm of 

spectroscopic techniques. 

The Sloan Digital Sky Survey (SDSS. lYork et al.ll200(tl 
presents a different model for studying elliptical galaxies. 
Unlike the surveys above that collect very detailed dynam- 
ical information for a small number of galaxies, the SDSS, 
with its 7r steradians of deep multicolour imaging and 10 6 
spectra, will measure global parameters for a very large 
number of galaxies, allowing one to statistically approach 
questions of galaxy structure. Indeed, the properties of el- 
liptical galaxies in the SDSS have a lready been studied i n 
great detail (Bernardi et al l2003al . l2003bl 12003d. l2003dh . 
This work constructed a sample of 9000 galaxies from early 
SDSS data, and exhaustively measured the various correla- 
tions between different observables. Principal results include 
a Faber- Jackson relation between luminosity L and velocity 
dispersion <x, a oc ^°' 25±0 012 ) anc j a fundamental plane re- 
lation between the effective radius R, the half-light surface 
brightness I, and a, R oc ^±0-05 j-o.75io.cn Thig sta _ 

tistical approach has been further extended using the weak 
lensing information from galaxy-galaxy lensing at large radii 
jMcKav et al.ll200ll) and models satisfying the dynamical 
constraints from velocity dispersions and weak len sing in an 
average sense have been constructed (|Seliak 20q2). 

This paper explores the dark matter content of ellipti- 
cal galaxies in the SDSS in more detail. Is there evidence 
for dark matter in ellipticals? How does it correlate with 
the luminosity and size of the galaxy? How does the dy- 
namical mass compare to the stellar mass of elliptical galax- 
ies? The simplest virial prediction combined with a constant 
mass to light ratio implies that L oc a 2 R\ however, the ob- 
served fundamental plan e (FP) of elliptic als shows the scal- 
ing L oc o- 1 - 98 ^ - 66 feernardi et al.ll2003cl) . There have been 
a number of proposals in the literature to explain this de- 
viation from the virial prediction, or "tilt". Implicit in the 
virial prediction is the assumption that elliptical galaxies 
form a homologous sequence, i.e. the dynamical structure 
of elliptical galaxies are self similar and related by simple 
scaling relations. Deviations from such a sequence would 
naturally manifest themselves as a tilt in the FP. However, 
base d on detailed dynamic al measurements of nearby ellipti- 
cals, iGerhard et all d200ll) argue that there is little evidence 
for any deviations from homology. A second, and preferred, 
proposal is to assume that M/L is not constant, but varies 
with luminosity. Such a variation could be caused either due 
to a metallicity (or dust)- luminosity correlation, or due to 
an increase in the dark matter fraction with luminosity. Is 
it possible to distinguish between these scenarios? These are 
the questions that this paper will attempt to answer. 

We start in Section 2 with the criteria for selecting our 
sample of SDSS elliptical galaxies. In Section 3, we propose 
a mass model for these galaxies and use it to estimate dy- 
namical masses based on measured scale sizes and velocity 
dispersions. We then compare stellar and dynamical mass 
estimates as a function of galaxy properties and discuss 
their implications (Sections 4 and 5). The appendix sum- 
marizes the relevant properties of the Sersic profiles we use. 
Wherever needed, this paper uses (fi m ,f2A) = (0.3,0.7) and 
H = 70 km/s/Mpc. 



2 THE SAMPLE 

The SDSS is imagin g 10 4 deg 2 of the N orthern Galactic 
Cap in 5 bandpasses iFukugita et alJll996|. u. q. r. i. z ) using 
a drift scanning, mosaic CC D camera jGunn et all Il998lf) 
under photometric conditions dHogg et "lu^|2"o"o" l|) . and is tar - 
geting 10 6 objects for spectroscopy (IBlanton et alj |2003b'l. 
most of which arc gala xies with r ba nd apparent magni- 
tude m r < 17.77 dStrauss et alJ 12002') . The data used in 
this paper cover an area of ~ 2000 deg 2 and include ~ 
160,000 galaxies with spectra, and are denoted samplelO 
within the SDSS collaboration. All of these data have been 
reduced by highly auto mated photometric and s pectroscopic 
reduction pipelines fsee IStoughton et al J l20o3 . for details). 
The astrometric calibration is automatically performed by 
a pipeli ne that obtains a bsolute positions to better than 0.1 
arcsec iPier et alJl2003h . and magnitudes are calibrated to 
a standard star ne twork approximately in the AB system 
femith et al]l2002D . 

The SDSS pipelines return a wealth of information for 
all detected objects; in addition, a number of auxiliary pa- 
rameters have been measured by various members of the 
SDSS collaboration for these objects. The principal param- 
eters relevant to this work are mentioned below. 

• Redshifts(z) and Velocity dispersions(cr): Each of the 
SDSS galaxy spectra is fit to a linear combination of galaxy 
templates at varying redshifts, broadened by a Gaussian ker- 
nel (Schlegel et al. I2003T) . Minimizing x 2 over this suite of 
models leads to an estimate of both the redshift and the 
stellar velocity dispersion of the galaxy. 

• Petrosian Fluxes : The primary measure of galaxy flux 
in the SDSS is the Pet rosian magnitude , a modification of 
the qu antity defined bv IPetrosianl dl976lk see lStrauss et alJ 
(2002) for details. Note that in the absence of seeing, the 
Petrosian flux is about 81.5% of the total flux of a deVau- 
couleurs profile. The photometric pipeline also returns the 
radii that enclose 50% and 90% of the Petrosian flux, R Pt 5o 
and Rp,go respectively, as well as the ratio of the minor to 
major axes, a/b, of the galaxy from model fits to a deVau- 
couleurs profile. 

• K-corrections : All luminosities and colours used in this 
paper are k-corrected to the median redshift of the survey, 
Zmed = 0.1, using the kcorrect package of IBlanton et al] 
j2003aD . 

• Sersic Profiles : The radial intensity profiles of all of the 
galaxies in our sample hav e been fit to Sersic profiles dSersid 
Il968t IBlanton et aljliool . 

i(R) = A 0xp [-(e/e a ) 1/n ] , (i) 

where R s is a scale radius, and n measures the concentration 
of the intensity profile . Unlike the Petrosian magnitudes, 
the Sersic fits are convolved with the PSF. Therefore, the 
effective or half light radius of the Sersic profile, -R50, is a 
robust estimator of the galaxy size and is what we adopt 
throughout this paper. 

• -P4000 : The c ontinuum break at 4000 A, D4000 
(Balo gh et al.l ll999'l is one of the most prominent features 
in galaxy spectra, and is caused by the accumulation of a 

1 We summarize the relevant properties of the Sersic profile in 
the Appendix. 
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Figure 1. Properties of the sample of elliptical galaxies. The 
upper left panel shows the distribution of the logarithm of stellar 
masses, while the upper right panel is the distribution of the r- 
band absolute magnitudes. The lower left panel shows the Sersic 
indices for our sample, where n = 4 is the de Vaucouleurs profile, 
and the lower right panel is the distribution of the g-r k-corrected 
colour. 



large number of metal lines in a narrow region of the spec- 
trum. Hot stars, indicative of a young stellar population, 
show a weak 4000 A break as the principal metals causing 
the absorption are multiply ionized, making this feature a 
powerful age es timator. The str ength of the break is com- 
puted following iTreniontll (|2002F ) (see also IKauffmann et all 

• Stellar masses (M*): IKauffmann et alJ (I2003T) compute 
stellar mass to light ratios for all spectroscopically observed 
galaxies by fitting stellar population synthesis models to 
the measured -D4000 st r ength and the HSa absorption line, 
assuming the iKroupal (1200 ll) initial mass function (IMF). 
Given this M/L ratio, the stellar mass is estimated by multi- 
plying by the Petrosian luminosity in the z band, correcting 
for extinction due to dust by matching the predicted galaxy 
colours to the observed colours. 



2.1 Selecting Ellipticals 

In order to translate the velocity dispersions into dynamical 
mass estimates, it is important that we select galaxies that 
are pressure supported, and not rotationally supported. In 
practice, this closely corresponds to selecting a sample of 
elliptical galaxies. We note that our selection criteria will 
also allow SO's; these galaxies, although not strictly ellipti- 
cal galaxies, satisfy our criterion that they have little or no 
rotational support. However, for simplicity, we refer to all 
these galaxies as elliptical galaxies. Since we desire minimal 
contamination from spirals, we use a combination of spec- 
tr oscopic and photometri c criteria, similar in spirit to those 
in lBernardi et all i2003alh although differing in detail. Note 



that unless specified, all the photometric quantities are mea- 
sured in the r band. The cuts made were 

• D4000 > 1.6 : As mentioned above, a large D4000 
corresponds to an older stellar population, normally as- 
sociated with ell i ptical s. This particular cut is based on 
IKauffmann et alJ (12003th 

• Rp,9o/R P ,50 > 2.6 JShimasaku et al.l l200ll) : Elliptical 
galaxies have intensity profiles that are more concentrated 
than spirals, allowing one to separate them on the basis of 
a concentration parameter, R p ,go/ Rp,so- 

• No emission lines : The stellar populations of ellipti- 
cals are generally old, with little or no star formation. Since 
emission lines are associated with star forming regions, we 
reject galaxies with emission lines (Ha, Nil) detected with 
S/N > 3. 

• a > 70 km/s : This cut is based on the resolution of 
the SDSS spectra; velocity dispersions lower than 70 km/s 
cannot be reliably determined. 

• Minor/Major axis ratio (a/6) > 0.7 : This was cho- 
sen to eliminate large edge-on spiral galaxies that were not 
eliminated by the previous cuts. Also, since a large rotation 
velocity would tend to flatten elliptical galaxies, this elimi- 
nates galaxies with large rotation velocities. 

Making these cuts reduces our sample from 165,812 to 
29,469 galaxies. Fig. shows some of the properties of our 
sample. Most of our galaxies have Sersic index n between 3 
and 4 and have their Petrosian g — r colour narrowly dis- 
tributed around ~ 1.0. We did not cut on either of these 
quantities, yet they lie in the range expected for elliptical 
galaxies. In addition, we have visually inspected images of 
a random subsample of these galaxies and have found the 
contamination from spirals to be less than 5 %. 



3 A DYNAMICAL MODEL 

In order to relate the velocity dispersions to dynamical 
masses, we need to develop a model of the mass distribu- 
tion of the galaxy. We model each galaxy (Fig|5| top panel) 
as having two spherically symmetric comp onents, a stellar 
component described by a lHernauisd <ll990T) profile, 



~\r + aY 



(2) 



and a dark matter halo. We adopt the Hernquist profile as it 
provides a convenient analytical approximation to the depro- 
jected de Vaucouleurs profile 2 . Modelling the dark matter is 
more difficult since it involves understanding its response to 
the b aryons. We start with an NFW profile jNavarro et alJ 
119971) , 



p oc r 



_1 (r + r 3 



(3) 



for the initial matter distribution, assume that a fraction F 
(the stellar mass fraction) condenses into stars, and that 
the remaining matter (the dark matter mass fraction) is 
adiabatically compressed by the stellar matter following 



2 See the Appendix for an explanation of why this is still justified 
for the Sersic profile fits that we use. 



© 0000 RAS, MNRAS 000, 000-000 



4 Padmanabhan et al 



iBlumenthal et all ||l986j). Determining the final dark mat- 
ter profile involves solving the equations 

nM NF w{n) = r f (M,(r f ) + M DM (r s )) 

(1 - F)M NFW {n) = M DM (r f ) , (4) 

where MNFw(fi) is the initial mass profile, M*(r/) and 
Mom{i"{) are the final stellar and dark matter profiles, and 
F is the stellar mass fraction. The first equation represents 
angular momentum conservation, while the second enforces 
the fact that orbits of dark matter particles do not cross. 
We now use the Je ans equation to find the 3 D velocity dis- 
persion profile, Vr l|Binnev fc Tremainelll987i) . 

r ydmr dmr J 

where M(r) is the total (stellar and dark) mass within r, v c 
(the circular velocity) is the velocity a particle on a circu- 
lar orbit would have, and /3 measures the anisotropy of the 
velocity distribution, 

/j=i-4 • ( 6 ) 

While the applicability of this model to elliptical galax- 
ies remains unclear, numerical simulations suggest it gives 
a good qualitative description in the region of interest 
(>0.5i?5o) (M. Steinmetz, private communication). Fig. [5] 
(top panel) shows the initial and final integrated dark mat- 
ter mass as a function of radius, as well as the stellar mass. 
The final circular velocity profile shown in Fig. [5] (solid line, 
bottom panel, the velocity was reduced by 1.65 as discussed 
below) is nearly flat over the range of interest, in agreemen t 
with observations of normal ellipticals jQerhard et all200ll) . 
The reader is referred to Seliak (2002) for more examples 
and the circular velocity profiles of individual components. 

We then can project this distribution to two dimensions, 

mw) -»jT(. -4) A- <7 » 

where I(R) is the 2D intensity distribution. The SDSS spec- 
trograph measures the luminosity weighted average (within 
the fiber aperture) of the above quantity, 

ff R' dR' I(R')a 2 (R') 

°( R ) = • (8) 

/„ R'dR'I{R') 

Examples of velocity dispersion profiles for different values 
of f3 are shown in the lower panel of Fig. [5] 

To compute the dynamical mass, we must choose a char- 
acteristic radius and relate a at that radius to the circular 
velocity, v c . The first complication is that the anisotropy 
j3 is unknown; however, as Fig. [5] suggests, the effects of (3 
are more pronounced in the inner regions suggesting that we 
measure a at large radii. This is a manifestation of the virial 
theorem; for an isothermal profile, the luminosity weighted 
velocity dispersion is simply related to the circular velocity 
via v\ = 3a 2 . From Fig.|2] we see that i> c (i?5o) ~ l-65a(Rso) 
with a scatter of ~ 10%, depending of the anisotropy profile 
we use. This gives us a dynamical mass estimate at R50, 

_ (LgScrfggO (Q) 

Mdynm = — ■ [V) 

Gr 




100 t 1 

0.5 1 1.5 2 2.5 3 

Figure 2. The upper panel shows the integrated mass of the 
original matter (dotted/red) profile, the stellar (solid/black) and 
final dark matter (dashed/blue) distributions. The lower panel is 
the luminosity weighted velocity dispersion profile for an isotropic 
= (solid), j3 = 0.3 (lower dotted) and = 0.5 (upper dotted) 
velocity distribution. We also show a variable anisotropy profile 
(dot-dashed), = 40ox/(l + x) 2 where x = r/rp, and /3q = 0.5 
and rp = 1.38 kpc. The short dashed (red) line is the circular 
velocity profile reduced by 1.65, while the long dashed (blue) line 
is the profile that would be observed by a slit spectrograph (using 
P = 0). Both panels assumed a galaxy halo of 10 12 A/q with 
a concentration (r v i r ial/ r s) of 10, and a stellar component of 
W il MQ with an effective radius of 5 kpc. 



This dynamical mass estimate is a 3D mass, while the stel- 
lar masses are projected masses. We estimate the 3D stellar 
mass by considering the mass contained within the Hern- 
quist profile at -R50, the projected half-light radius, which 
is approximately ~ 42% of the total stellar mass. Further- 
more, the stellar masses are computed using the Petrosian 
flux, which is approximately 80% of the flux of an early type 
galaxy; these two factors together imply that the 3D stellar 
mass at Rz,q is ~ 50% of the total estimated stellar mass. 

Unfortunately, the velocity dispersion estimated by the 
SDSS spectra is not at -R50, but at the fiber diameter of 3 
arcseconds. One solution is to follow IrBernardi et alJ (|2003a) 
and apply an empirical correction to the measured disper- 
sion. Another approach, also suggested there, is to stack 
galaxies of similar masses and physical sizes and construct 
a velocity dispersion profile. These galaxies are at a range 
of redshifts, so the 3 arcsecond fibers probe the velocity dis- 
persion at different physical radii, allowing us to construct 
a composite velocity profile. We can then use this profile to 
correct the measured velocity dispersions out to Rho- 

However, naively making subsamples is dangerous be- 
cause the SDSS spectroscopic sample is magnitude limited, 
i.e. to a good approximation, it included all galaxies with 
m r < 17.77. This implies that galaxies at higher redshifts are 
more luminous on average, which could introduce spurious 
correlations between observed quantities. Therefore, to en- 
sure uniformity, all the subsamples we construct in this pa- 
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Figure 3. Estimated velocity profiles for our sample of ellipticals 
in SDSS as a function of Rf/R^o where Rf is the fiber radius 
of 1.5 arcseconds. The different panels are different stellar mass 
volume limited subsamples, where the panels are labelled by the 
logarithm of the stellar mass range. The different line segments 
in each panel show the median velocity dispersion of subsamples 
in the physical size of the galaxy. Note that the size of the galaxy 
increases from right to left. The various size ranges considered are 
(from top to bottom) 2.5-6.0, 3.0-6.0, 4.0-10.0 and 5.0-13.0 kpc. 
The volume limited subsamples ensure a uniform sample over the 
redshift range of interest. The velocity profiles derived here are 
consistent with being flat over the region of interest. 



per are volume limited, including only those galaxies which 
would remain in the sample at the redshift limits. 

Fig. shows the velocity profiles estimated by stack- 
ing galaxies. The different panels are different volume lim- 
ited mass samples that uniformly sample the redshift range, 
while the different line segments are the median velocity 
dispersions of subsamples of similar physical sizes. The pro- 
files are consistent with being flat. This would appear incon- 
sistent with the velocity dispersion profiles measured from 
indi vidual ellipticals w hich, in general, decline with radius 
(e.g. iKronawitter et al.ll2000t) . However, the profiles are sig- 
nificantly flatter when luminosity weighted within the fiber, 
shown in Fig|5| comparing the profile that would be seen 
through a slit with that through a fiber. We therefore use 
the velocity dispersion measured through the 3 arcsecond 
aperture as the dispersion at i?so. While the statistical ve- 
locity dispersion profiles are in a good agreement with the 
isotropic model predictions, the narrow dynamic range of 
volume limited subsamples prevents us from making any 
strong conclusions on the anisotropy parameter f3. 



4 RESULTS 

Now that we have both stellar and dynamical masses, we can 
compare them. Figs. [I] and |K| show the stellar and dynamical 
mass to light ratios respectively as a function of luminosity. 
A useful consistency test of both the stellar and dynamical 
masses is that Mdynm/L (~ 5 — 8) is greater than M*/L 



Figure 4. The upper panel is the stellar mass to light ratio 
(M t /L) as a function of r band luminosity. The different seg- 
ments/colours correspond (from left to right) to volume limited 
subsamples within redshift intervals 0.04-0.08 (black), 0.08-0.12 
(red), 0.12-0.16 (blue), 0.16-0.20 (green), 0.20-0.24 (magenta). 
The solid line is the median of the distribution, while the lower 
and upper dashed lines are the 16% and 84% intervals. The lower 
panel shows the distribution of M-*/L at the median luminosity 
of each redshift subsample, staggered by A log(M* /L) = 1 for 
clarity. 




Figure 5. Same as Fig. HI except that the dynamical mass to 
light (M dynm / L) and the dynamical mass to stellar mass ratio 
(M t j !/nm /M*) are plotted. Although we don't show it explicitly 
here, the scatter is consistent with log-normal. 
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(~ 3). Since the stellar masses are sensitive to the choice 
of IMF, this constrains the possible IMF's as we discuss in 
the next section. Also, M dynm /L increases roughly as L 017 , 
while M*/L is constant, or slightly decreasing. The apparent 
decrease in M*/L with luminosity could be interpreted as 
luminosity evolution, although the statistical significance of 
the decrease is unclear (see below). Note that the scatter in 
both quantities is consistent with being log-normal, with a 
width of ~ 30% for the stellar masses and slightly more for 
the dynamical masses. 

A relevant question is whether this scatter is due to 
measurement error or an intrinsic scatter in galaxy proper- 
ties. The errors in the stellar masses are obtained by inte- 
grating over a grid of population synthesis models to obtain 
the 16% and 84% confidence intervals; these errors are ~ 
25%. Estimating an error for the dynamical mass is more 
involved; we assume a 10% error on the conversion of a to 
the circular velocity (see Fig|5]for a justification) and a fur- 
ther ~ 10% error on the measurement of a 2 Rso- This yields 
an error in the dynamical mass of ~ 30 %. This implies that 
the dominant source of scatter is measurement errors, so it is 
possible that the intrinsic correlations are tighter than what 
is observed. 

The measurement errors can be dominated by system- 
atic or statistical errors. Systematic errors can be caused by 
the velocity dispersion-mass conversion we have assumed, er- 
rors in extrapolation to R50, or by errors in the stellar mass 
determination. Pure statistical errors, such as those from 
measurement errors, will average out in median quantities. 
If statistical errors dominate, then the trends in the medi- 
ans of these distributions are highly statistically significant. 
If the errors are dominated by systematics, the trends may 
still be significant, as long as the systematic error does not 
couple to the parameter that is varied. In all our figures, 
we have presented the average Iff contours for individual 
objects. This is the most conservative error estimate and it 
should be kept in mind that the errors on median quantities 
may be significantly smaller than this. 

However, there are possible uncertainties in the stel- 
lar mass estimations that could lead to systematic errors. 
The stellar masses are estimated by comparing the mea- 
sured D4000 and HSa indices with detailed stellar population 
synthesis models. In order to estimate the robustness of the 
masses to these details, we co mpared the stellar m asses used 
here with those computed bv lPanter et al.l <|2002T) and found 
that, for the subsample of ellipticals found in both catalogs, 
the differences were consistent with the errors. Specifically, 
all the trends found in thi s paper are a l so pre sent for the 
stellar masses estimated bv lPanter et alJ (|2002). 

Another potential source of systematic error is colour 
gradients: M*/L is determined from the spectra within the 
fiber, typically at 0.4-0.8i?5o, and then assumed to hold 
within i?5o. We can try to measure an average M„/L pro- 
file using the same method we used to estimate the veloc- 
ity dispersion profile. This is shown in Fig|S| as before, the 
different panels are different stellar mass ranges, while the 
different line segments are the median M*/L for galaxies 
of similar physical sizes. There is a modest trend towards 
decreasing M*/L with radius, although the statistical sig- 
nificance is unclear. If we attribute the observed trends to a 
colour gradient, our stellar masses are only upper limits on 
the true stellar mass. The suggested trends would reduce the 
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Figure 6. M*/L profiles for our sample, analogous to Fig|3 The 
mass subsamples are labelled in each panel, while the individual 
lines are medians for different physical sizes. 

stellar mass by at most 20%. While this would alter the pre- 
cise form of the correlations that we have observed, it would 
increase the dark matter fraction observed, but would not 
change any of the trends we find. However, luminosity evo- 
lution would also cause a similar trend and is difficult to 
disentangle from M*/L gradients. Therefore, we choose to 
be conservative and assume that the estimated M*/L value 
is applicable at R$o- As a final check on systematics, we find 
that our trends remain unchanged even if we select only 
very spherical objects with a/b > 0.9 that are unlikely to 
show any rotation. This provides additional support to our 
assumption of ignoring rotation in the analysis. 

We can now consider galaxies of a given luminosity and 
ask if there is any variation in Md ynm j 'M* as a function of 
the size of the galaxy. In order to do this, we start with the 
volume limited samples discussed above and construct lumi- 
nosity subsamples 0.5 magnitudes in width, centred about 
the median luminosity of each original sample. Fig.|7|shows 
Mdynm/M, as a function of R50, for samples created in this 
manner. The ratio of dynamical mass to stellar mass in- 
creases roughly linearly with radius, with a slope of ~ 0.14. 
Fig. |H| shows the same relation for stellar and dynamical 
mass subsamples; the relations for these subsamples are very 
similar to those for luminosity. 

It is interesting to ask whether these correlations can 
be explained with the model of Section [3] The model has 
four parameters, the virial mass and scale size of the halo, 
and the stellar mass and scale size of the galaxy. The latter 
two of these are directly constrained by observations. We 
parametrise the scale size of the halo (Eq|3J by the concen- 
tration parameter, c = r V i r /r s , where the virial radius, r„ir> 
is determined by the virial mass, 

47T -i 

M mr = —Apcritrtir , (10) 

where A(= 200) is the spherical overdensity and p cr u is the 
critical density. While numerical simulations suggest that 
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Figure 7. Mdynm/ 'M* as a function of the physical size (half 
light radius) of the galaxy. The different segments/colours cor- 
respond to cutting the sample in redshift as in Fig. [3] and then 
selecting a narrow luminosity bin around the median luminosity 
of the redshift subsample. The median luminosities for the sub- 
samples are (from left to right) are -19.93, -20.78, -21.30, -21.78 
and -22.20. 



the concentration decreases with virial mass, this variation 
is small over the range considered here. We thus make the 
assumption that the concentration is constant (c = 12.7) 
independent of the size of the halo. This leaves us with only 
the virial mass of the halo to vary; Fig.[5]compares the pre- 
diction of this model for the dynamical to stellar mass ratio 
as a function of R 5 q with the observed result, for galax- 
ies with stellar masses ~ 1.5 x 10 n M Q . The model does 
well predicting the general trend for this and other stellar 
masses, although the slopes are shallower than what is ob- 
served. While the statistical significance of the disrepancy 
is small, it is possible that it is caused by the restrictive as- 
sumptions of the model. Allowing either the concentration 
or virial mass to vary with the size of the galaxy of a given 
stellar mass can reproduce the observed slope. 

This simple one-parameter family of models provides a 
novel method of estimating the virial mass of the halo in 
which a galaxy of a given stellar mass resides. This method 
is demonstrated in Fig|5] one adjusts the virial mass to best 
match the observed correlation of Mdynm /M, with R50, 
while "errors" can be estimated by matching to the 16% 
and 84% contours. The results of this exercise are in Fig llUI 
as a function of stellar mass. Here again we take the con- 
servative approach of assuming the errors to be dominated 
by systematics. If the errors are statistical, then the usual 
1/vN factor can be applied and the errors on the mean re- 
lation become considerably smaller. This figure argues that 
the virial masses of halos are 10 times the stellar mass of 
ellipticals at the high mass end and somewhat less at the 
lower mass end. 



Figure 8. Same as FiglTI except that the subsample is cut in 
stellar (upper panel) and dynamical (lower panel) mass. From 
left to right, the median stellar masses are 4.36, 9.33, 15.3, 22.8 
and 31.8 x10 10 Mq, while the median dynamical masses are 3.57, 
9.51, 16.5, 27.1, 42.2 xl0 10 M Q . 

5 DISCUSSION 

In this paper, we have compared stellar and dynamical mass 
estimates for a sample of elliptical galaxies in the SDSS, and 
studied the dependence of that relation on the luminosity 
and size of the galaxy. 

The first conclusion that we draw from this work is that 
stellar mass estimates are consistent with the dynamical es- 
timates, i.e. M* < Mdynm with A/* ~ Mdynm for high sur- 
face brightness (low R50) galaxies. This is far from a trivial 
conclusion, since the procedure to estimate stellar masses 
is a rather involved one. In particular, the stellar masses 
are extremely sensitive to the stellar init i al ma ss function. 
Her e we have follo wed iKauffmann et aD i2005tl and adopt 
the iKroupal feOOll) IMF; replacing it with a Salpeter IMF 
increases the number of low-mass stars, effectively increas- 
ing all the inferred stellar masses. An examination of Fig. 
Q shows that the stellar masses would no longer be consis- 
tent with the dynamical estimates in that case for high sur- 
face brightness (low R50) galaxies. This statement depends 
on how one regularises the divergence of the Salpeter IMF 
at low masses, and therefore, must only be interpreted as a 
constraint on the minimum mass that the Salpeter slope can 
extend to. Our results suggest that a cutoff at M > O.IMq 
is required, consist ent with IMF measu rements from young 
embedded clusters jMuench et alj |2002l. 

In what follows, we assume that the difference between 
the dynamical and stellar masses is due to the presence of 
a dark matter component of the galaxies. Ellipticals often 
have a hot gas compone nt; however, it is a relatively minor 
contribution to the mass llSarazi rJl997ft . and the stellar mass 
is, to a good approximation, the baryonic mass of the galaxy. 

Given that approximation, we can draw the following 
conclusions about elliptical galaxies in SDSS : 

• The fraction of dark matter increases with increasing 
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Figure 9. The solid and dotted lines show the median and 16% 
and 84% contours of the distribution of Mdynm/M* as a function 
of R50 for a volume limited subsample with stellar mass between 
1.17 and 1.86 xlO 11 M (median value of 1.5 X 10 n M Q ). The 
central dashed line shows the prediction of the model in Sec 131 
for the median stellar mass and a halo mass of 1.35 X 1O 12 M , 
while the upper and lower lines have halo masses of 6.0 X 10 11 Mq 
(lower) and 4.8 X 10 12 Mq. All the initial halos had concentration 
parameters c = 12.7. The vertical dotted line is the median size 
of all the galaxies in the sample. 



a 




5xl0 10 



Figure 10. The ratio of the virial mass to stellar mass as a 
function of stellar mass calculated using the methods described 
in the text and Fig[Tj] The shaded region is excluded by the cosmic 
baryonic fraction of ~ 17 % measured by WMAP, assuming only 
baryons within the virial radius can cool and form stars. 



luminosity: As is evident from Figs0] and |KJ the stellar 
mass to light ratio is approximately independent of luminos- 
ity, while Mdynm/L oc L 017 , implying an increasing dark 
matter fraction with luminosity. Note also that the depen- 
dence of Mdynm/L on L is consistent with the results of 
iBernardi et all (l2003bl) . M dynm /L oc L 014±0 02 . 

This result suggests that the tilt in the fundamental plane 
is caused by an increasing dark matter fraction with lumi- 
nosity. Recall that the virial prediction for the fundamental 
plane, 



(tHw 



oc RI , 



(11) 



is traditionally modified by assuming that the mass to light 
ratio varies with luminosity, [M/Li] oc L a . As discussed in 
the introduction, there are two possibilities to explain this 
scaling, variations in metallicity or in the dark matter frac- 
tion. Variations in metallicity would show up as variations 
in Mt/L, while variations in the dark matter fraction would 
be evident in Mdynm/L. We observe that M*/L is approxi- 
mately independent of L, while Mdynm/L increases with L, 
favouring the conclusion that the tilt is due to dark matter 
and not m etallicity. This is the op posite conclusion from that 
reached bv lGerhard et al] fcOOll) who argue that a maximal 
stellar mass is supported by population synthesis models, 
and therefore, one need not invoke dark matter. Our ap- 
proach is different; we use stellar masses estimated indepen- 
dently of the dynamical estimates. Comparing these shows 
that elliptical galaxies have a significant dark matter com- 
ponent within an effective radius, with as much as 3-4 times 
the stellar component for the largest galaxies. 

We can al s o co mpare our results with those of 
iBorriello et al] i2003T) . which attempted to fit the funda- 
mental plane with stellar and dark matter mod els similar 
to those we use. However, IBorriello et al] i2005t) leave the 
stellar mass to light ratio a free parameter; this leads them 
to conclude that the data favours a model with no dark 
matter and M*/L ~ 5.3. This is approximately a factor of 
2 greater than our estimates of M*/L ~ 3; moreover, we ar- 
gue that changing our IMF to agree with the higher M*/L 
would make the stellar masses inconsistent with our dynam- 
ical estimates for high surface brightness galaxies. 

• The dark matter fraction increases with increasing size 
at constant luminosity (Fig. [7J . Note that this is a robust 
result, since similar trends are also seen if the stellar mass or 
velocity dispersion are kept constant, instead of the luminos- 
ity (Fig|HJl. For the highest surface brightness galaxies, the 
stellar mass approaches the dynamical mass. This suggests 
that the IMF based on local observations is also applicable 
to ellipticals. A simple model of a galaxy embedded in a dark 
matter halo with a constant star formation efficiency 3 cap- 
tures the qualitative trends. In such a model a low surface 
brightness galaxy and a high surface brightness galaxy of 
equal luminosity sit in equal dark matter halos, but the low 
surface brightness galaxy is more spread out, so the ratio of 
stellar to dark matter mass is lower. These trends are thus 
expected in ab-initio models of galaxy formation. To explain 
them without the dark matter would require introducing an 
IMF that depends on the surface brightness. 



3 In other words, that the stellar mass fraction, F, is independent 
of the size of the galaxy. 
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While this model can explain most, if not all, of the trends 
suggested by the observed relations, one question remains: 
what sets the scale size of the galaxy? As Fig. |7| shows, the 
half light radius varies by a factor of two at any given lu- 
minosity. One simple possibility is that it is related to the 
scale radius of dark m atter, which simulati ons suggest has 
a considerable scatter feullock et al.ll200ll) . However, if a 
more compact galaxy is in a more concentrated dark matter 
halo then this would make the agreement in Fig. |5] worse, 
since that figure suggests we need less dark matter in com- 
pact galaxies and more in low surface brightness galaxies. 
In this case the assumption of star formation efficiency be- 
ing independent of galaxy scale size must also be violated. 
While we do not have definitive answers to these questions, 
we note that ongoing galaxy-galaxy lensing studies on the 
same SDSS sample may be able to provide useful additional 
information, since it can determine the virial mass as a func- 
tion of luminosity, surface brightness etc. 
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We also use this model to estimate the virial mass as 
a function of mass and luminosity for our sample of galax- 
ies. Although the errorbars are large, we find that the virial 
mass of the halo is around 10 times the stellar mass of the 
galaxy (Fig. I1UI . We note that our modelling assumes that 
halos are undisturbed. Many of the ellipticals are in denser 
enviroments of groups and clusters. If the galaxy is at the 
center of the group then our modelling still applies. If it is a 
satellite then some of the dark matter attached to the galaxy 
before it merged into the group halo was likely stripped off 
due to the tidal effects. In this case our estimates apply to 
the virial mass of the galaxy prior to merging. This is the 
mass of interest if most of the stellar mass has been assem- 
bled prior to the galaxy merging into the group or cluster. 

Comparing 10% stellar mass frac tion to the cosmic 
baryonic fraction of 17% from WMAP JSpergel et alJl2003l) 
implies that 60% of the initial gas has been converted to 
stars; a la uncertainty takes us between 25% and 100%. 
This is not necessarily inconsistent with the global stel- 
lar mass fraction, which is only around 5-10% of baryons 
jBaloeh et alJl200ll: iDickinson et al.ll2003l) : it suggests that 
the efficiency of star formation peaks at M v i T ~ 10 12 Mq and 
is likely to be much smaller in either lower and higher mass 
halos. We observe this trend at the high mass end, where 
the dark matter fraction increases with increasing mass, al- 
though the uncertainties in our measurement are large. The 
implication of an increasing dark matter fraction with lumi- 
nosity is that the star formation is suppressed in the most 
massive halos, as expected for example by the longer cooling 
times in hotter halos. The sample of low mass ellipticals is 
too small to observe the reverse trend at the low mass end. 
Both the absolute value of the star formation efficiency and 
its trend with luminosity /stellar mass are in agreement with 
galaxy-galaxy lensing results, which provide an independen t 
estimate of the total mass of the halo (iGuzik fc Seliakl2002D . 
These additional constraints on the mass also allow for a 
more detailed modelling of the structure of ellipticals. Ul- 
timately, we hope that these results will both inform and 
constrain theories of galaxy formation. 

We would like to thank Daniel Eisenstein, David Hogg, 
Raul Jimenez, Yeong-Shang Loh and Scott Tremaine for use- 
ful discussions. U.S. acknowledges support from the Packard 



APPENDIX A: 
PROFILES 



PROPERTIES OF SERSIC 



The Sersic profile is an obvious generalisation of the classi- 
fication of galaxies into deVaucouleurs and exponential in- 
tensity profiles, 



I(R) oc exp(— x 



l/ns 



(Al) 



where a; is a dimensionless radial variable, R/R s - Exponen- 
tial and deVaucouleurs profiles correspond to n = 1 and 
n = 4 respectively; in general, the Sersic index n can be 
any positive real number. Below, we summarise the relevant 
properties o f Sersic profi les; a more detailed discussion is 
contained in Ciotti (1991) and references therein. 

The first quantity of interest is the half light radius R50 . 

we can solve for the integrated intensity 



Defining y \ 
profile, 

L(R) «/» dy'y' 2 "- 1 exp(- 
oc 7(2n,y) , 



(A2) 



where 7(m, y) is the incomplete Gamma function. Using the 
above result, we see that R50 can be determined by solving, 



7(2n,ft(n)) 
T(2n) 



(A3) 



where F(2n) = 7(2n, c») is the Gamma function and b(n) = 
{R50/R3) 1 ■ This equation can be solved numerically, but 
a useful approximation is b(n) — 2n — 0.324. 

We can also solve for the deprojected Sersic profile v(r); 
this is obtained via the Abel transform, 



. , 1 r di 
v (r) = -- ^5 



dR 



dR 



(A4) 



In order to quantify the differences between the different 
values of n, it is useful to consider the total mass/light con- 
tained within a radius r. This is shown in Fig lAll where we 
have scaled all the radii by -R50. Note that the discrepan- 
cies between different values of n are marginal, especially at 
radii greater than 0.3i?so (which contains less than 20% of 



© 0000 RAS, MNRAS 000, 000-000 



10 Padmanabhan et al 




0.1 1 10 

r/R M 



Figure Al. Integrated light profiles for deprojected Sersic pro- 
files for n=3 [dotted, red], 4 [solid, black], 5 [short-dashed, blue], 
6 [long-dashed, green] profiles, as a function of r/R^o where iJso 
is the projected half-light radius. The vertical dashed line marks 
the 3D half light radius for the n=4 profile. 

the total mass/light). This justifies our modelling the stel- 
lar light distribution of all galaxies with a Hernquist profile, 
even when n/4. 
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